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Abstract 

It is shown that the threshold productions of n^pp, vr+nj? and n^d from 
pp collisions can be consistently described by a model consisting of a pion 
s-wave rescattering and NN pair-terms of heavy-meson exchanges. The large 
difference between a^°^{pp —> and a^°^{pp —> vr^np) is understood from 

the orthogonality of the deuteron and the np scattering wave functions. In a 
calculation using the Paris potential to account for the initial and final NN 
interactions, it is found that the data can be best reproduced by using a soft 
ttNN form factor with = 650 MeV for a monopole form. This is consistent 
with an earlier study of pion production in the A excitation region. 
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With the development of the Cooler facility at the Indiana University Cyclotron Facility 
(lUCF), accurate data of threshold pion production from proton-proton (pp) scattering have 
become available. The data of pp ppn^ []^ had led to the identihcation of the short-range 
pion production mechanisms |^|-^]. In this paper we demonstrate that by also including the 
new data of the pp —> n^np and pp —>■ reactions in a theoretical analysis, we can 

extract information concerning the range of the ttNN form factor. Furthermore, we will 
show that the large difference between the cross sections of pp —> and pp ir'^np can 

be understood from the behavior of the final np wave functions. 

Following the previous investigations M,0i, the s-matrix for the production of a pion 
with a charge a from a pp collision can be written as (we use the normalization of 0) 

Sfi = -2'Ki5^^\pi+p2-pf - ( 1 ) 

where is the initial pp scattering wave funtion, and (/| is either the deuteron wave 

function in pp n~^d reaction or the np scattering wave function {xp p\ in pp ^ ir^np 
reaction. Following the previous approach ||^, the plane-wave matrix element of the pro¬ 
duction operator is dehned by taking a nonrelativistic limit of the lowest-order Feynman 
amplitudes calculated from a model interaction Hamiltonian. For our present purpose, it is 
sufficient to consider the following interaction Hamiltonians for tt, u, a, and N fields 

H(x) = H^(x) + H,{x) + H^{x) (2) 

where 

f - 

H^{x) = — f'ip{x) ■ d^(j)^{x) 

Hs{x) = - - i){x)(j)^{x) ■ (j)^{x)ij{x) H-^ Vi(x)r V’(x) • 0(x) x —0 (x) 

\ ot 

Ha{x) = gai’ix)^p{x)(j)^{x) 

H^{x) = guj‘ip{x)'j^‘ip{x)df,(j)^{x) ( 3 ) 

Here describes the familiar pseudovector ttNN coupling. The vriV s-wave interaction is 
described by Hg, as introduced in Ref. [0. In this work, we use the vaules Ai ~ 0.005 and 
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\2 ~ 0.05 calculated p| from the ttN s-wave scattering lengths of |T^. The heavy meson- 
exchanges are induced by for a isoscalar-scalar a meson and for a isoscalar-vector uj 
meson. The isovector mesons 6 and p can be included , but will be neglected in this work 
for simplicity. 

The pion s-wave rescattering term of the production operator is calculated from 
H.^ and Hg. The NN pair-terms of a and uj exchange amplitudes are calculated from 
//tt and second orders in and H^. To ensure that there is no double counting in the 
calculation of Eq. (1), only the negative energy component of the nucleon propagaor in the 
heavy-meson Feynman amplitudes is kept. The resulting terms are called the NN pair- 
terms. To be consistent with the employed nonrelativistic NN potentials for calculating 
the initial and final NN wave functions, we make the usual nonrelativistic expansion of the 
Feynman amplitudes. In terms of momentum variables and pion isospin (iV(pi) -|- N{p 2 ) — 
^(Pi) + ^{p' 2 ) + '^{^1 ''^6 obtain the following form 


An, — 


+ (1 <—^2) 


( 4 ) 


where 


ri"=Gob(P2-p>2)(27r)! 


a {1) -k 




rp(TT) ^ q(-k) \ ^2 • {P2 — P' 2 ) 


—q‘^ — ml 




2m N 

(-l)“r_„(2) +i 


a(l) ■ (pi+j/i) 


, . 871 X 2 Uk + UJg 


m^r 2m, 


(-1)“ [r(i) X mu 


TU = <! ^ 2 "^" 2 \^U ■ {Pi+^1 + 


.q. _m-^ myr 


rr{uj) _ 

-‘-a ^ 

with 




- 9 I 


-g 2 L 


a(l) ■ (P2 +^' 2 ) + ij{p'2-P2) ■ ( 0 ^( 1 ) X <j(2)) 


(5) 


n — i 1 / 

0 {2npl'^ mn ’ 


( 6 ) 


Here we set q = P 2 - P 2 , ^k = \Jml + kX , ujg = ^Jml + q^, r± = ^ (n ± i T 2 ), tq = r^, 
and a is the z-component of pion isospin. The form factors F^'^\Ai,U) are introduced to 
regularize meson-NN vertices. 
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Let us first consider the pion s-wave rescattering term Its main feature is that 
Ai ~ 0.005 -C A 2 — 0.05. Since the matrix element of the isospin operator (ri x 72 ) between 
two T = 1 states vanishes, only the very weak Ai term of can have a contribution to the 
pp —>■ ppn^ cross section. That is why the calculations including only failed 

to account for the total cross section of pp ppn^. It was found that the discrepancies 
with the data can be removed by introducing the NN pair terms of heavy-meson exchange. 
A more general approach was introduced in Ref. [@] to calculate the NN pair-terms from 
realistic NN potentials. For the present purpose, it is sufficient to use the simple model 
dehned in Eqs. (5). We choose the monopole form factor ^. All of the 

two nucleon wave functions in Eq. (1) are calculated from the Paris potential [0. We found 
that the threshold pp —> ppir^ total cross section data can be best reproduced if we choose 
gl/An = 6 and = 1000 MeV/c, for the cr-exchange, and gll4:TT = 10 and A^ = 1400 MeV 
for the cn-exchange. In considering only the total cross sections, one can also ht the data by 
keeping only the a-exchange term with g'^/^i^ = 12 and A^- = 2000 MeV/c. This is similar 
to the Endings of Refs. and [Q. All results presented below are from the cx-exchange 
model. The results for a -I- cu-exchange model are very similar and are therefore omitted. 

Once the heavy meson-exchange mechanism has been determined in the study of the 
pp —> pp-K^ reaction, the only freedom of the model is the form factor FA)(^At^^ g^) of 
Since the T = 0 final np states are involved in the charged pion production, the large A 2 term 
of will contribute to the cross section. Its contribution depends on the range A^^ of the 
irNN form factor. We find that all of the threshold pion data can be best reproduced with 
A^ = 650 MeV/c. Our results are the solid curves displayed in Fig. 1. In the same figure, 
we also show the results (dotted curves) without including the NN pair-terms of heavy 
meson-exchanges. The heavy meson-exchange mechanisms clearly dominate the pp —>■ ppir^ 
cross section; but also plays a significant role in charged pion productions. 

An interesting feature of Fig. 1 is that at energies near the n~^np production threshold 
(292.30 MeV), the cross sections of the tt+c? production are almost two orders of magnitude 
larger than the vr+np production. This can be understood as follows. In this near threshold 
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energy region, the final np pair is mainly in the ^So{T = 1) and — ^Di{T = 0) states. 
The production of the ^So{T = l)np state is very weak mainly because the large A 2 term of 
the pion rescattering does not contribute (as in the case of pp —> ppn^). The difference 
between the 7T~^d and vr+np productions is therefore only in their final radial wave functions of 
the np system. With the same ^Si — ^Di quantum number, the np scattering wave function 
should be orthogonal to the deuteron wave function. As illustrated in Fig. 2, this comes 
about from any realistic NN potential. The orthogonality is due to the change of the sign 
of the np scattering wave function in the region where the deuteron wave function is large. 
Consequently, the vr+np production amplitude involves a cancellation in the integration over 
the production operator, and hence is very much suppressed. The results shown in Fig. 1 
reflect the correctness of the employed Paris potential in describing this wave function effect. 

We have found that our predictions are very sensitive to the range A^r of the form factor 
FC). This is not surprising since the threshold pion production involves a large momentum 
transfer at the nNN vertex. In Fig. 3, we show that the calculations with A^ = 1000 MeV/c 
overestimate significantly the charged pion production cross sections. This is consistent with 
the results from the earlier work [|^| in the A-excitation region. 

The results presented above are from the calculations using the ttNN coupling constant 
P/4:7r = 0.079 of Ref. [T^. In recent ttN |T^ and NN phase shift analyses, a smaller 

value of P/4:n = 0.075 was obtained. We have also carried out calculations using this 
smaller nNN coupling constant. A slightly larger ttNN cutoff A = 680 MeV is needed to 
obtain an equally good fit to all of the data shown in Figs.l and 3. 

In conclusion, we have shown that all of the threshold pion production data from pp 
collisions can be described by a model consisting of the pion s-wave rescattering term and 
the NN pair terms of heavy meson-exchanges. The large differences between the 7T~^d 
and 7r’''?7,p production are shown to be due to the orthogonality between the np scattering 
wave function and the deuteron bound state. Because of the isospin character of the pion 
rescattering term, the pp —>■ ppvr® data can be used to pin down the NN pair terms of the 
heavy meson-exchanges. The charged pion production can therefore be used to determine 
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the 71NN form factor. A soft ttNN form factor with = 650 MeV for a monopole form 
is found to be consistent with the data. To end, we mention that the present investigation 
is based on the distorted-wave formulation, Eq.(l), which was used in all of the previous 
studies |@-^,|^,|[. To make further progress, it is necessary to consider the coupling between 
the hnal vr+np and vr+d channels. This is a highly nontrivial task. Our effort in this 

direction will be reported elsewhere. 

I would like to thank Steve Dytman and Ralph Segel for informing me of their data, 
Vincent Stoks and Iraj Afnan for useful discussions. This work is supported by the U.S. 
Department of Energy, Nuclear Physics Division, under contract W-31-109-ENG-38. 
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FIGURES 

FIG. 1. The total cross sections (solid curves) calculated from using the Paris potential and 
= 16, Ag- = 2000 MeV/c for the cr-exchange model of are compared with the data. The 
dotted curves do not include the NN pair-terms of heavy meson exchanges. 

FIG. 2. The solid curve is the ^Si deuteron bound state wave function. The dotted and dashed 
curves are the ^Si scattering wave function for = 1 and 10 MeV. The wave functions are 
generated from the Paris potential. 

FIG. 3. The solid (dotted) curves are the calculated total cross section using the cutoff 
Ayr = 650 (1000) MeV for the nNN form factor. The results are from calculations using the 
Paris potential and = 16, = 2000 MeV/c for the cj-exchange model of Aa- 
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